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Abstract 

Numerical modelling of coherent spin relaxation in nanomagnets, formed by mag- 
netic molecules of high spins, is accomplished. Such a coherent spin dynamics can be 
realized in the presence of a resonant electric circuit coupled to the magnet. Com- 
puter simulations for a system of a large number of interacting spins is an efficient tool 
for studying the microscopic properties of such systems. Coherent spin relaxation is 
an ultrafast process, with the relaxation time that can be an order shorter than the 
transverse spin dephasing time. The influence of different system parameters on the 
relaxation process is analysed. The role of the sample geometry on the spin relaxation 
is investigated. 
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1 Introduction 



A polarized magnet placed in an external magnetic field, whose direction is opposite to 
the sample magnetization, comprises a strongly nonequilibrium system. The relaxation of 
spins to their equilibrium position can occur in different ways. The simplest process is the 
slow relaxation during the time Ti caused by the spin-lattice coupling [1], when the total 
magnetization tends to zero. This is an incoherent process, since Ti is usually much longer 
than the spin dephasing time T2. A straightforward way to make the spin relaxation coherent 
would be by imposing a strong transverse field simultaneously driving all spins, which would 
represent the process of free induction [1], lasting for the time T2. This, however, is a rather 
trivial process, when spins are practically independent of each other. 

A more elaborate situation arises if the magnet is inserted into a magnetic coil of a 
resonant electric circuit. Then the magnetic field, induced by the coil, provides an efficient 
feedback mechanism organizing the coherent motion of spins [2]. It is possible to realize 
six different regimes of coherent spin relaxation in a polarized sample coupled to a resonant 
circuit: Collective induction, maser generation, pure superradiance, triggered superradiance, 
pulsing superradiance, and punctuated superradiance. A detailed description of these regimes 
is given in the review articles [3-5]. 

Collective induction and maser generation with spins were realized in a number of ex- 
periments [6-9]. Pure spin superradiance was, first, observed in Dubna [10,11] and later 
confirmed by the groups in St. Petersburg [12,13] and Bonn [14]. Pulsing superradiance was 
demonstrated by experiments in Ziirich [15,16]. The regime of punctuated superradiance 
was suggested in Ref. [17] and, to our knowledge, has not yet been realized experimentally. 
A comprehensive survey of experiments can be found in Refs. [4,5,18]. 

It is necessary to emphasize the principal role of the resonant electric circuit, coupled to 
the spin sample, for realizing the regimes of spin superradiance. This makes the fundamental 
difference between the spin superradiance and the atomic superradiance [19]. The latter 
can be achieved in a resonatorless system [19-21], though a resonant cavity can enhance the 
effect [22,23]. Contrary to this, the spin superradiance, occurring in the radiofrequcncy range, 
cannot be realized without a resonator, which is due to the destructive role of the dipole spin 
interactions and to the absence of the feedback mechanism collectivizing the spin motion. 
This basic difference was emphasized in Ref. [19] and thoroughly explained in Refs. [5,24,25]. 
One should also distinguish between coherent transient effects, caused by intense alternating 
external fields, and having much in common for optical atomic systems [20,21], gamma 
radiation [26-29], as well as for spin samples [30-32], and superradiant phenomena, when 
the self-organization of coherence is the basic origin of the arising superradiance [25,33,34]. 

Here we concentrate our attention on the superradiant regime of spin motion, which re- 
quires the presence of the resonant electric circuit, providing the feedback mechanism for the 
collective self-organization of spin motion. The regimes of free induction, collective induction, 
and triggered spin superradiance could be considered in the frame of the phenomenologi- 
cal Bloch equations supplemented by the Kirchhoff equation for the circuit [35-38]. These 
phenomenological equations, however, are not applicable for describing pure spin superra- 
diance, when the coherent motion of spins develops in a self-organized way from initially 
chaotic spin fluctuations. The complete theory, based on the microscopic spin Hamiltonian, 
and describing all regimes of spin relaxation has been developed in Refs. [18,25,33,34,39] 
and expounded in detail in the review articles [4,5]. This theory is in good agreement with 
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experiment [4,5] as well as with numerical computer simulations [40-42] accomplished for 
nuclear or electron spins one half, S = 1/2. 

There also exists a wide class of paramagnetic materials formed by magnetic molecules 
[43], whose effective ground-state spins can reach rather high values of 5* ~ 10. These 
molecules compose molecular magnets whose properties are described at length in Refs. 
[5,18,44,45]. For example, the molecules Mni2 and Feg possess the spin S = 10. The 
magnetic cluster compound Mn604Br4(Et2dbm)6 has the total spin S* = 12 [46]. And the 
magnetic molecule Cr(CNMnL)6(C104)9, where L stands for a neutral pentadentate ligand, 
displays the effective spin S = 27/2 [45]. 

Molecular magnets, formed by high-spin magnetic molecules, exhibit a rather strong mag- 
netic anisotropy and a gigantic relaxation time of their magnetization, which reaches about 
two months in zero magnetic field at low temperature [5,18,44,45]. The magnetization relax- 
ation is mainly due to the phonon-assistcd mechanism caused by spin-phonon interactions 
[47,48]. The fast superradiant-type relaxation of magnetization is theoretically possible for 
molecular magnets, provided they are placed in a sufficiently strong external magnetic field 
and necessarily coupled to a resonant circuit [5,18,25,49]. 

The main points, distinguishing the relaxation in magnetic molecules, studied earlier 
[44,45,50-57], from the process to be investigated in the present paper, are as follows. 

First of all, the usually studied spin relaxation in magnetic molecules is the process of 
quantum tunneling occurring in individual molecules [50-57]. Contrary to this, we shall be 
interested in coherent, principally collective effects, due to strong correlations between many 
molecules, which is caused by the resonator feedback field. 

Quantum tunneling in separate molecules is known to be phonon assisted, with an essen- 
tial phonon infiuence on the relaxation process [47,48,58]. But in the case of spins strongly 
correlated through the resonant feedback field, spin-phonon interactions are not as impor- 
tant, the basic dynamics being governed by the resonator feedback field. 

The single-molecule spin relaxation is not connected with a fixed particular frequency 
[59], because of a sufficiently strong magnetic anisotropy and a varying external magnetic 
field [5,19,24,25]. While to realize the collective coherent spin relaxation necessary requires 
the situation close to resonance between the Zeeman frequency and the resonator natural 
frequency. 

Magnetic avalanches are usually accomplished by sweeping the longitudinal magnetic 
field with a rather slow sweeping rate. A typical field sweeping rate used in experiments 
[50-55] is about 0.1 T/s. In our case, we consider a fixed external magnetic field, with a well 
defined Zeeman frequency. 

The avalanche time of the magnetic moment in magnetic molecules is between 10~^ s 
and 1 s [50-57,59] . This is rather slow process, as compared to the coherent spin relaxation, 
which is an ultrafast process, with characteristic relaxation times about 10~^° — 10~^^ s 
[5,18,24,25,49]. 

The peculiarities of the fast superradiant-type spin relaxation in molecular magnets can 
be successfully analyzed by means of computer modelling, which, to our knowledge, has 
not yet been accomplished. This approach provides a very efficient tool for studying the 
microscopic properties of spin system. And it is the aim of the present paper to describe the 
results of computer simulations for analyzing the features of the fast coherent relaxation of 
molecular spins, typical of the high-spin molecular magnets. 

The outline of the paper is as follows. In Sec. II, we present the main definitions and 
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equations to be employed in our computer modelling. In Sec. Ill, the results for a bulk 
sample are analyzed. Since the dipolar spin interactions are anisotropic, it is interesting to 
study the related anisotropic geometric effects for different shapes of the magnetic samples. 
These geometric effects are investigated for a chain of molecular spins oriented either along 
the external magnetic field or perpendicular to it and for spin planes, with the external 
magnetic field being either perpendicular to it or lying in that plane. Section IV contains 
conclusions. 



2 Basic Definitions and Equations 

We consider a spin sample characterized by the molecular spin vectors Sj = {Sj,Sj,Sj} 
associated with the lattice sites enumerated by the index j — 1,2,. . . ,N. An external 
magnetic field is directed along the z-axis. 

Bo = Boe,. (1) 

This defines the Zeeman frequency 

<^o = - ^ -Bo = ^ )Ub -Bo , (2) 

in which /xq = — is the electron magnetic moment, with /is being the Bohr magneton. In 
general, similarly to the magnetic-resonance setup, there can exist the transverse magnetic 
field, directed along the x-axis, 

Bi = Bie^ , Bi = ho + hi cos cut , (3) 

and consisting of a constant field ho and an alternating field hi cos ut. 

Molecular magnets possess the single-site magnetic anisotropy characterized by the anisotropy 
parameter D, which defines the anisotropy frequency 

u;n = (2S-l)j, (4) 

where S is the molecular spin. The magnetic anisotropy exists for high spins, playing an 
important role, while for 5* = 1/2 it disappears, according to Eq. (4). 

Spins interact with each other through the dipolar forces characterized by the dipolar 
tensor 



< - ^, (^«/^ - <.) , (5) 

in which 



II ^ij 

fij = I ^ij I ) ^ij = ) r = Tj Tj . 

fij 



For what follows, it is convenient to introduce the dipolar coefficients 

a,,^D-, b,j^l{D^f-Dyy-2tD^y) c,j^l{D^f-zDyj) , (6) 
having dimensions of energy. 
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The spin sample is inserted into a magnetic coil of an electric circuit characterized by 
the circuit damping 7 and the circuit natural frequency uj. Moving spins generate electric 
current in the coil, which, in turn, produces the feedback magnetic field H acting on the 
spins. The generated electric current is described by the Kirchhoff equation. Choosing the 
coil axis along the x-axis, so that H = He^, the Kirchhoff equation can be rewritten [33,34] 
as the equation for the feedback magnetic field H, 

— + 2jH + I H{t') dt' = -47r77 , (7) 

where the effective electromotive force in the right-hand side of Eq. (7) is produced by 
moving spins, with the average magnetization 

// ^ 

- y E < 5, > , (8) 

V being the sample volume. The fiUing factor 77 in the right-hand side of Eq. (7) is approx- 
imately equal to = V/Vc, where is the coil volume. For what follows, without the loss 
of generality, we may assume the dense filling, with rj — 1. Instead of Eq. (7), we can use 
the equivalent differential equation 

(PH dH ^ d'^rrix 

_ + 2,-+.^H = -4.^. (9) 

in which we set 77 = 1. 

All possible attenuation mechanisms have been carefully described in Ref. [25]. Those 
that influence the spin motion are as follows. The longitudinal attenuation 

ri = 7i + 7i* (10) 

is the sum of the spin-lattice attenuation 71, caused by spin-phonon interactions, and of 
the polarization pump rate '-fl, due to a stationary nonresonant pump, if any. The total 
transverse attenuation is 

r2 = 72(1-5') +72*- (11) 

This includes the homogeneous broadening 72, renormalized by the factor (1 — s^), appearing 
in the case of strongly polarized spin systems [1,25], with s being the average spin polarization 
reduced to the number of spins N and to the spin value S. The last term 72 is the static 
inhomogeneous broadening. 

In order to represent the equations of spin motion in a compact form, it is convenient to 
introduce the ladder spin components 

SJ- = S] - iS^ , S+ = S] + iS^ . (12) 

Also, we shall use the following notation: 
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Ci^l i: hc.,S] - I <^^JS- + 2k,Sl) . (13) 
The effective force, acting on tfie j-spin, can be written as 

/,^-^/xo(5i + i/) + 0- (14) 

The derivation of equations of motion for the spin variables S*", S'^ ^ and 5*1 has been 
described in great detail in Refs. [4,5,25]. The resulting equation for reads as 

-^ = -1 (o^o + e° - ^r2) ^7 + j^s] + i ^ s]s^ . (15) 

The equation for 5*^ is conjugate to Eq. (15). And the equation for is 

lf = -5(/;s7 + s;/,)-r,(sj-c), (16) 

where is the stationary spin polarization. From Eqs. (15) and (16), with notation (12), 
one can always return to the evolution equations for 5^, ^J, and Sy 

In numerical simulations, one treats the spins Sj as classical vectors [40-42]. It is conve- 
nient to work with the reduced quantities characterizing the reduced transverse spin variable 

1 ^ 

and the reduced longitudinal spin variable 

The spin variables (17) and (18) characterize the collective properties of a large number of 
magnetic molecules composing the molecular magnet. The time evolution of these variables 
is prescribed by the equations of motion (15) and (16). This picture of collective spin motion 
is a generahzation of the evolution equations for a single magnetic molecule. The study of 
collective coherent effects is the main aim of the present paper. 

In our numerical simulations, we solve the spin evolution equations (15) and (16) for a 
finite number of spins N . The resonator feedback field is given by Eq. (9), with the initial 
conditions 

//(O) = , ij(0) = , (19) 

where the overdot imphes the time derivative of H. The spin variables at the initial 
time are distributed randomly over the sample, so that to obtain a prescribed value s(0) of 
the spin polarization (18), while variable (17), for sufficiently high initial spin polarization, 
being negligible, 

s(0) = So , ^^(0) = . (20) 

The external magnetic field (1), with Bq > 0, is aligned with e^. For the initial spin 
polarization sq > 0, the magnetic moment of the molecular sample 

M(0) = NSfioSoe, = -2NSfiBSoe, 

is opposed to e^. That is, the considered molecular magnet is prepared in a strongly nonequi- 
librium initial state, from which it relaxes to a stationary state. 
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3 Results of Computer Simulation 



First, we consider the coherent spin relaxation in bulk samples, where spins are located 
in lattice sites of a cubic lattice. Computer simulations are accomphshed for N — 125 
spins. For larger values of N, the results are quahtatively similar. The periodic boundary 
conditions have been imposed. Since our aim is to study the self-organized process, we 
set zero the transverse field Bi = 0. We assume that the spin sample is without defects, 
so that the inhomogeneous broadening is negligible, 72 = 0. For low temperatures, the 
spin-lattice interaction in molecular magnets is very weak, with the longitudinal relaxation 
time Ti reaching months (see review articles [5,18,44,45]), because of which the attenuation 
parameter Fi plays no role, and we can set Fi = 0. It is convenient to deal with dimensionless 
quantities measuring all frequencies in units of 72, so that we set 72 = 1. And we shall 
measure time in units of 7^^, that is, in units of T2. 

First and foremost, we have to stress the necessity of the resonant circuit. When the 
latter is absent, there is no fast relaxation at all. Then there could exist only the very 
slow polarization decay during the time T2, which is caused by spin interactions. The 
same slow relaxation happens if there is a circuit, but there is no resonance between its 
natural frequency lu and the Zeeman frequency ujq. Therefore, in what follows, we always set 
(J — coq. This resonance condition is necessary, though not sufficient. To reahze an effective 
spin reversal, it is important that cuq be much larger than the anisotropy frequency (4). The 
condition ujq ujd ensures that the anisotropy does not induce an effective detuning from 
the resonance [5,18,24,25]. 

Figure 1 illustrates how the spin reversal depends on the value of the Zeeman frequency 
ujq. The larger the latter, the more pronounced is the spin reversal. 

The resonator damping 7 defines the resonator ringing time r = I/7, during which the 
magnetic sample effectively interacts with the resonator. The relation between the resonator 
ringing time r and the transverse relaxation time T2 essentially infiuences the spin reversal. 
When r = T2, then there is a permanent exchange of energy between the spin sample and 
the resonator, so that the spin polarization oscillates around zero. When r = 0.1 T2, there 
occurs a well pronounced reversal, hence the value 7 = 10 is optimal for the latter. And 
if T = 0.02 T2, then the effective interaction time between the sample and resonator is 
too short to realize a good reversal of polarization. The corresponding three qualitatively 
different cases ate shown in Fig. 2. 

The magnitude of spin reversal also depends on the initial polarization. The larger the 
initial value so, the stronger the spin reversal, which is illustrated in Fig. 3. 

Magnetic anisotropy is an obstacle for the coherent spin relaxation. The larger the value 
of ujo, the smaller the spin reversal, as is demonstrated in Fig. 4. 

Dipole spin interactions is also a factor suppressing spin coherence. This is illustrated 
by Fig. 5, where the behavior of spin polarization for the case with dipole interactions 
is compared with that one, for which the dipole interactions are switched off by setting 
D"^!^ = 

It is interesting that switching off the dipole interactions yields the figures that are very 
close to those obtained by the reduction of spin from 5* = 10 to 5* = 1/2. Thus the dashed 
line in Fig. 5, where 7 = 10, for 5" = 10 can be repeated not by setting the dipole tensor 
to zero but by reducing the spin to 5" = 1/2. In Fig. 6, we show the behavior of spin 
polarization for 5" = 10 and 5" = 1/2 for 7 = 1. Again, switching off the dipole interactions 
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for 5" = 10 yields the dashed curve corresponding to S = 1/2 with dipole interactions. 

Dipole interactions are anisotropic. Therefore, their influence on the relaxation process 
can be different for the samples of different shapes and of different orientations with respect 
to the external magnetic field and with respect to the resonator feedback field. To analyze 
these differences, we study the spin relaxation, under the same system parameters, but for 
different samples. We consider the chain of spins oriented either along the external magnetic 
field, i.e., along the 2;-axis, or along the feedback field, that is, along the x-axis. And we 
consider the plane of spins, oriented either in the y — z plane or in the x — y plane. The results 
of computer simulation for = 144 spins are presented in Fig. 7. These results demonstrate 
that the maximally efficient spin reversal happens for the chain of spins directed along the 
X-axis. This looks quite understandable, since the x-axis is the axis of the direction of the 
resonator feedback field, which is the main source of the coherent spin motion. 

4 Discussion 

We have accomplished computer simulations of the coherent spin relaxation in molecular 
magnets with large spin. The investigation is based on the microscopic model taking into 
account realistic dipole spin interactions and the single-site magnetic anisotropy. The system 
is prepared in a strongly nonequilibrium state, with an external magnetic field opposite to 
the sample magnetization. 

The principal point of our investigation is the presence of a resonator coupled to the 
sample. The later is inserted into a coil of an electric circuit, whose natural frequency is in 
resonance with the Zeeman frequency. Without the resonator, the coherent spin motion is 
impossible. It is the resonator feedback field, which collectivizes the spin motion, making it 
well correlated, hence, coherent. 

To realize the coherent spin relaxation, the Zeeman frequency uuo has to be much larger 
than the anisotropy frequency ljd. An efficient spin reversal requires that the initial spin 
polarization be sufficiently high, the higher, the better. The typical spin reversal time r is 
an order smaller than the transverse dephasing time T2, which translates into the relation 
7 ~ IO72. 

The role of dipole spin interactions, in the presence of a resonator, is twofold, making 
the spin dynamics in a sample coupled to a resonator rather different from that happening 
in a sample with no resonator feedback fields. 

Prom one side, dipole interactions infiuence the spin motion by making the spin reversal 
less pronounced. For low spins, such as 5 = 1/2, dipole interactions are less important than 
for large spins S = 10. Emphasizing the decoherence infiuence of the dipolar interactions, 
we should keep in mind that our simulations are performed for a finite number of spins. The 
majority of our calculations are done for 125 spins. Because of the long range of the dipolar 
forces, increasing the number of spins strengthens the decohering infiuence of these forces. 
However, all qualitative results remain, as we have checked by varying the number of spins 
between 64 and 343. Also, presenting the results in dimensionless units, as we have done, 
when all frequencies and attenuation parameters are normalized by the dipolar interaction 
strength, makes the calculated curves for the average magnetization practically independent 
of the sample size. Increasing the number of interacting spins simply implies the renormal- 
ization of the dimensionless quantities and does not change their behavior represented in 
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dimensionless units. 

At the same time, increasing the number of spins strengthens the role of the resonator 
feedback field, which makes the process of spin coherentization faster, so that the relaxation 
time, due to the coherent spin motion, depends on the number of spins as 1/N. 

In this way, stronger dipole interactions, from one side, increase the transverse decoher- 
ence attenuation, but, from another side, they induce stronger coherence through the action 
of the resonator feedback field, making the coherent relaxation faster. These two opposite 
effects, to some extent, compensate each other. Therefore, the coherent spin dynamics, oc- 
curring in the presence of a resonator, qualitatively does not change much under the variation 
of spin number. 

Our main concern in the present paper has been the study of spin dynamics for large 
spins. This is why we have done numerical simulations for S = 10. We have had no aim of 
studying the low-spin dynamics, such as that of spins one-half, since this dynamics has been 
considered earlier. It is only to note that the low-spin dynamics is rather different that we 
show its qualitative difference in one curve of Fig. 6. 

It is worth mentioning that for large spins jS" > 1 it is feasible to consider spin transitions 
between different sublevels labelled by the 2;-projection number m = —S, —S+1, . . . , S—1, S. 
A pair of sublevels can be treated as an effective two-level system [56,57]. Then to realize the 
coherent spin relaxation, one has to tune the resonant natural frequency to the transition 
frequency between the chosen two levels. For high nuclear spins, this procedure was realized 
experimentally [15,16]. Hence, in the same way it can be realized for molecular magnets. 

It is important to stress that there are several principle physical differences between the 
experiments with nuclear spins, described in Refs. [15,16], and the situation considered in 
the present paper. In these experiments [15,16], the nuclei of ^^Al inside the ruby crystal were 
studied. First of all, the nuclei of ^^Al possess spins 1 — 5/2, which are not as large as we 
have considered here, dealing with S = 10. Second, in the case of ^^Al, an external resonant 
circuit was tuned to the central line { — 1/2, 1/2}, with a fixed transition frequency uq ~ 10^ 
Hz, thus, reducing the consideration to an effective two-level system with spin one-half, while 
here we always have dealt with the total spin S = 10, since we have considered the resonant 
circuit tuned to the transition between —S and S. Third, as we have shown, for our high-spin 
case the influence of dipolar interactions is essential, while their role for an effective one- 
half spin system [15,16] is not of such importance. Fourth, contrary to the case of nuclear 
spins, having no single-site anisotropy, the molecular magnets, we have studied, exhibit 
quite strong magnetic anisotropy, fundamentally distorting spin dynamics and making it 
principally different as compared with the isotropic case of nuclear spins. Fifth, for strongly 
polarized spin materials, it is necessary to take account of the saturation effect making the 
total transverse attenuation depending on the polarization level [1,25], as in Eq. (11), while 
this effect does not play role for not so strong polarization [15,16]. Sixth, in experiments 
[15,16], pulsing spin dynamics was analyzed, when the inversion of spin polarization was 
permanently supported by constantly applied dynamic nuclear polarization with a rather 
high pumping rate, while we studied the pure coherent spin relaxation, when there is no 
permanent pumping. Finally, we have studied here the dependence of spin relaxation on 
the sample shape and orientation. Such geometric effects, to our knowledge, have not been 
investigated. These seven factors make the spin dynamics in our case and in the case of 
Refs. [15,16] principally different. 

When studying the geometric effects related to the sample shape and its orientation. 
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differently oriented spin diains and planes have been considered. We have found that, 
under the same system parameters, including the number of spins, except the sample shape 
and orientation, the most efficient coherent spin relaxation, with the deepest spin reversal, 
happens for a chain of spins aligned with the direction of the resonator feedback field. 

The main aim of the present paper has been to analyze the coherent spin relaxation under 
widely varying system parameters, in order to clarify the influence of different parameters on 
the coherent spin motion. This should help in choosing the optimal materials for realizing 
such a coherent spin motion. Nowadays, there are plenty of magnetic materials, with widely 
varying properties, which could be used for experimentally observing the described effects. 
The description of the properties of various molecular magnets can be found in the review 
articles [5,18,44,45]. 

As an example, we can mention the most often studied molecular magnets made of the 
molecules Mni2 or Fcg, whose spins are S = 10. For these materials, uju ~ 10"*^^ s^^. At low 
temperatures, below the blocking temperature of about 1 K, the sample can be polarized, 
having a very long spin- lattice relaxation time Ti ~ 10^ — 10^ s. Hence 71 = 1/Ti is 
practically neghgible, 71 ~ 10~^ — 10~^ s~^. Dipole spin interactions in these materials are 
rather strong, with 72 ~ lO^*' s~^. To realize the coherent spin relaxation, the external 
magnetic fleld Bq is to be sufficiently strong, such that corresponding Zeeman frequency ujq, 
being close to the resonator natural frequency tu, would be much larger than the anisotropy 
frequency (Vd- For the considered case of Mni2 or Fes, this requires the field Bq ~ 100 T. 
This is a strong field, though which can experimentally be reached [60]. Fortunately, there 
are many other molecular magnets with a smaller anisotropy. For instance, in the case of 
nanomagnets formed by the molecules Mug, whose spin is S* = 12, the magnetic anisotropy 
is much lower, with ljd ~ 10^" s~^. Therefore the required external magnetic field is only 
about So ~ 1 T, which is the standard field used in laboratory. 

The existence of the magnetic anisotropy typical of many molecular magnets, hinders 
the feasibility of the coherent spin relaxation. However, there are many magnetic materials 
with a small anisotropy, which should not disturb the coherent spin motion. In addition, 
the influence of the magnetic anisotropy can always be suppressed by a sufficiently strong 
external magnetic field. Fortunately, there exists in the world the possibility of creating 
very strong magnetic fields. Among available sources [60] , we may mention those where the 
magnetic fields up to 45 T (USA) and even 600 T (Japan) can be reached. 

In order to estimate the typical time of the coherent spin relaxation, we may notice that 
this time is an order smaller than T2. The spin dcphasing time in molecular magnets, such 
as Mni2 and Fes, is due to dipole interactions yielding 72 ~ lO^'^ s~^. Hence T2 ~ 10~^° s. 
This means that the typical time of the coherent spin relaxation in these materials is 10~^^ 
s, which is an ultrafast process. 

In conclusion, it is worthwhile to mention that for many magnetic molecules the influence 
of hyperfine interactions from nuclear spins could be important. These interactions result 
in the appearance of an additional line broadening, which can be inchidcd in the effective 
attenuation parameters, so that the existence of the hyperfine interactions can be taken 
into account by the appropriate definition of the effective attenuations, as has been done in 
Ref. [42]. At the same time, for many molecules, typical of the large family of magnetic 
molecules, such as Mni2 and Fcg, the hyperfine interactions are of the order of 10~^ K, which 
are much weaker than the dipolar interactions, being of the order of 0.1 K (see details in 
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the review articles [5,18,44,45]). When the hyperfine interactions are two orders smaller 
than the dipolar interactions, the former, with a very good approximation, can be neglected. 
And if the former arc comparable with the latter, this can be taken into account by the 
corresponding definition of the line broadening. 
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Figure Captions 



Fig. 1. Reduced spin variable s, for a cubic lattice, characterizing the spin polarization 
along the z-axis, as a function of dimensionless time (measured in units of T2) for the 
Zeeman frequencies coo — 1000 (solid hne), loq — 2000 (long-dashed hne), and loq — 5000 
(short-dashed line). The simulation is done for the molecules of spin S — 10, with the 
reduced initial polarization sq — 0.9. The anisotropy frequency is cud — 20 and the resonator 
damping is 7 = 10. 

Fig. 2. Reduced spin polarization s, as a function of dimensionless time, for a cubic 
lattice, with uq = 2000 and ujd = 20, for the resonator damping 7=1 (solid line), 7 = 10 
(long-dashed line), and 7 = 50 (short-dashed line). The sample of molecules with spin 
5" = 10 has the initial polarization sq = 0.9. 

Fig. 3 Reduced spin polarization s, as a function of dimensionless time, for a cubic 
lattice, with uq = 2000, ud = 20, 7 = 10, S = 10, for different initial polarizations sq = 0.9 
(solid line), Sq — 0.7 (long-dashed line), and Sq — 0.5 (short-dashed line). 

Fig. 4 Spin polarization s, as a function of dimensionless time, for a cubic lattice, with 
loq — 2000, 7 = 10, 5" = 10, and for different magnetic anisotropy values characterized 
by the anisotropy frequency cud — 20 (solid line), cud — 50 (long-dashed line), cud — 100 
(short-dashed line). 

Fig. 5. Spin polarization s, as a function of dimensionless time, for a cubic lattice, 
with loq — 2000, lod = 20, 7 = 10, and S — 10, for two different cases, when the dipole 
interactions are present (solid line) and when they are absent (dashed line). 

Fig. 6. Spin polarization s, as a function of dimensionless time, for a cubic lattice, with 
loq = 2000, cuz) = 20, 7 = 1, and for different spins 5" = 10 (sohd hne) and S — 1/2 (dashed 
line) . 

Fig. 7. Difference in the behavior of spin relaxation for different sample shapes and 
orientations, under the same values cuq = 2000, ujd = 20, 7 = 30, S = 10. The chain of spins 
along the z-axis (solid line) ; the chain of spins along the x-axis (long-dashed line) ; the plain 
of spins in the y — z plane (short-dashed line), and the plane of spins in the x — y plane 
(dashed-dotted line). 
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Figure 1: Reduced spin variable s, for a cubic lattice, characterizing the spin polarization 
along the z-axis, as a function of dimensionless time (measured in units of T2) for the 
Zeeman frequencies ujq = 1000 (solid line), ujq = 2000 (long-dashed line), and Uq = 5000 
(short-dashed line). The simulation is done for the molecules of spin S = 10, with the 
reduced initial polarization Sq — 0.9. The anisotropy frequency is lud — 20 and the resonator 
damping is 7 = 10. 
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Figure 2: Reduced spin polarization s, as a function of dimensionless time, for a cubic lattice, 
with ujo — 2000 and lod — 20, for the resonator damping 7=1 (solid line), 7 = 10 (long- 
dashed line), and 7 = 50 (short-dashed line). The sample of molecules with spin 5" = 10 has 
the initial polarization sq — 0.9. 
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Figure 3: Reduced spin polarization s, as a function of dimensionless time, for a cubic lattice, 
with Uq = 2000, uj£, = 20, 7 = 10, S* = 10, for different initial polarizations Sq = 0.9 (solid 
line), So = 0.7 (long-dashed line), and Sq = 0.5 (short-dashed line). 
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Figure 4: Spin polarization s, as a function of dimensionless time, for a cubic lattice, with 
loq — 2000, 7 = 10, 5" = 10, and for different magnetic anisotropy values characterized 
by the anisotropy frequency cud — 20 (solid line), cud — 50 (long-dashed line), cud — 100 
(short-dashed line). 
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Figure 5: Spin polarization s, as a function of dimensionless time, for a cubic lattice, with 
ujq = 2000, U£) = 20, 7 = 10, and S = 10, for two different cases, when the dipole interactions 
are present (solid line) and when they are absent (dashed line). 
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Figure 6: Spin polarization s, as a function of dimensionless time, for a cubic lattice, with 
uq = 2000, ujd = 20, 7 = 1, and for different spins 5* = 10 (solid line) and S = 1/2 (dashed 
line) . 
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Figure 7: Difference in tfie beliavior of spin relaxation for different sample shapes and ori- 
entations, under the same values uoq = 2000, = 20, 7 = 30, S = 10. The chain of spins 
along the 2;-axis (solid line); the chain of spins along the x-axis (long-dashed line); the plain 
of spins in the y — z plane (short-dashed line), and the plane of spins in the x — y plane 
(dashed-dotted line). 
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